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We investigate the initiation and subsequent evolution of Crow instability in an inhomogeneous 
unitary Fermi gas using zero-temperature Galilei-invariant non-linear Schrodinger equation. Con- 
sidering a cigar-shaped unitary Fermi gas, we generate the vortex-antivortex pair either by phase- 
imprinting or by moving a Gaussian obstacle potential. We observe that the Crow instability in a 
unitary Fermi gas leads to the decay of the vortex-antivortex pair into multiple vortex rings and 
ultimately into sound waves. 
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I. INTRODUCTION 



A vortex-antivortex pair in a three dimensional invis- 
cid fluid is susceptible to the long wavelength coopera- 
tive instability known as Crow instability (CI) [1 . In 
an undisturbed vortex-antivortex pair, each vortex feels 
the velocity field produced by the other. As a result of 
it, the pair moves in the direction perpendicular to the 
plane containing the pair. Now, some initial disturbance 
in the fluid can lead to the formation of the sinusoidal 
distortions in the vortex shape, which are symmetric with 
respect to the plane separating the two vortices. In this 
case, the total velocity fleld experienced by the vortex 
(antivortex) has three contributions: (a) velocity fleld 
induced by the vortex (antivortex) on itself due to its 
curvature, (b) the straining experienced by the vortex 
(antivortex) when displaced from its equilibrium position 
in the velocity field induced by the undisturbed antivor- 
tex (vortex), and (c) the velocity field induced on the 
vortex (antivortex) by the sinusoidally distorted antivor- 
tex (vortex). The strain rate is directly proportional to 
circulation of the vortex and inversely proportional to 
square of the distance between the vortices. The first 
and the sum total of the last two contributions are also 
termed as the velocity fleld produced due to the self and 
mutual-induction respectively. In the absence of the sec- 
ond vortex and hence mutual-induction, the sinusoidal 
displacement disturbance travels along the vortex line 
and rotates about its unperturbed position, thus distort- 
ing the vortex line into a helical shape. These helical 
waves, which are neutrally stable, propagating along the 
length of the vortex line are also known as Kelvin waves 
[2]. In the presence of the second vortex, the self in- 
duced rotation can be decreased or increased due to the 
azimuthal velocity component of strain fleld. When the 
two balance each other, the unbalanced radial component 
of the strain fleld leads to the exponential growth of the 
displacement perturbation; this is known as CI. As the 
amplitude of these sinusoidal oscillations reaches a criti- 
cal value, the vortex and antivortex reconnect leading to 
the formation of a chain of vortex rings. It was shown by 
Crow that a vortex line in the presence of only strain fleld 
will be unconditionally unstable, whereas it will be neu- 



trally stable of only self induction were present [T]. The 
interplay between these two opposing effects leads to the 
CI. The CI is also responsible for the decay of wing-tip 
generated vortex-antivortex pair in the wake of the air- 
craft T]. The CI in this case can be triggered by the 
turbulence in surrounding air or local variation in the air 
temperature or density. It was shown by Kuznetsov et al. 
[3], that only the long- wavelength symmetric modes are 
unstable due to CI. The long-wavelength antisymmetric 
modes, on the other hand, are stable [3]. The typical 
wavelength of the unstable modes is much larger than 
the size of the vortex core [T]. 

The CI can occur not only in classical fluids but also 
in superfluids. The characteristic feature of the superflu- 
ids is the quantized circulation in them 01 Ej- Hence, the 
realization of quantized vortices in Bose- Einstein conden- 
sates (BECs) [SH^ and two-component Fermi superfluids 
[S| has established the superfluid nature of these systems. 

In case of homogeneous BECs, the CI has been in- 
vestigated theoretically in Refs. [31 HO]- In this context, 
Kuznetsov et al., [3] demonstrated the long wavelength 
instability of acoustic solitons and vortex pairs in the 
frame work of three dimensional Schrodinger equation. 
It was noted in this study that in case of solitons, the 
instability corresponds to snake instability [TT], while in 
case of vortex-antivortex pairs it takes the form of CI. Re- 
cently, CI in cigar shaped BECs was theoretically studied 
in Ref. [T^]. On the experimental front, the snake insta- 
bility of the solitons leading to the formation of vortex 
rings has been observed [T31 [Tl] , whereas the CI is still 
eluding the experimental observation. The CI is inhib- 
ited in oblate BECs. This is due to the suppression of 
bending and reconnections of vortices in oblate BECs 
[15] . Also, the generation of Kelvin waves jUlIT^, which 
aids the development of the CI, is suppressed in oblate 
condensates [TH]. 

In case of two-component Fermi gas, the attractive in- 
teraction between the two component spin states is a 
prerequisite for superfluidity. Depending upon the value 
of s-wave scattering length and Fermi wave vector fcp, 
these superfluids can be weakly interacting /cplal ^ 1 
or strongly interacting fcpre ^ 1 < feplaj, where 
is the effective range of interatomic interaction poten- 



2 



tial. The two limiting cases, kp\a\ — )• (a — > —0) and 
fcplal — >■ CJO (a — )■ — oo) are known as Bardeen-Cooper- 
SchriefFer (BCS) and unitary limits respectively. The 
inverse of Fermi wave vector fcp is of the order of aver- 
age interparticle separation. Thus, the effective range is 
much less than the interparticle separation which in turn 
is much less than s-wave scattering length in the unitary 
limit. It implies that the unitary Fermi gas (UFG) is si- 
multaneously dilute (ro ^ fcp^) and strongly interacting 
{kp^ ^ \a\). In the unitary Fermi gas (UFG), the ther- 
modynamic quantities become the universal functions of 
Fermi energy Ep and the ratio of temperature to Fermi 
temperature T/Tp [TH]. The generation, real time evo- 
lution, and interaction among the quantized vortices has 
been studied in UFG using the time-dependent superfluid 
local density approximation at zero-temperature |20j . 
The method involves the computationally intensive task 
of solving time dependent coupled Bogoliubov-de Gennes 
equations self-consistently. An alternative method based 
on the single orbital time dependent density functional 
theory has been developed in Refs. [?TI - I55] to study the 
superfluid Fermi gases. Employing this approach, the 
UFG is well described by a non-linear Schrodinger equa- 
tion (NLSE) at T = OK, which leads to same superfluid 
density as the original many body system. The NLSE 
is appropriate to study UFG, provided the characteris- 
tic wavelength of the phenomenon under study is larger 
then the coherence length [5S] . The NLSE has been used 
to study various phenomena like collective modes p]"H23] , 
expansion dynamics [HJI^llllS], solitons [151 HZ]) Joseph- 
son oscillations [55], shock wave generation [551 [57]. etc. 
Recently, NLSE has also been used to investigate the soli- 
ton mediated decay of shock waves into vortex dipolcs in 
the oblate BCS superfluids 38J. In addition to all these, 
the Bose-Fermi superfluid mixtures can be described by 
the Gross-Pitaevskii equation (GPE) for the bosons, cou- 
pled to the NLSE for the order parameter of the Fermi 
superfluid [57| [SD]- Such coupled equations have been 
used to study gap solitons [57], mixing-demixing tran- 
sition , the spontaneous symmetry breaking of the 
Bose-Fermi mixture in double- well potentials [40], local- 
ization of a Bose-Fermi mixture in a bichromatic optical 
lattice [Hj , etc. The NLSE with quintic repulsive nonlin- 
earity has also been derived to describe one dimensional 
Tonks- Girardeau (TG) gas [42], and was used to study 
the formation of shocks and their dynamics [13]. The 
same equation, with an additional nonlocal cubic term 
accounting for the dipole-dipole attraction, has been used 
to study existence, stability and dynamics of bright soli- 
tons in dipolar one dimensional TG gas [44 . 

In the present work, we theoretically study the CI in 
strongly interacting Fermi superfluids using NLSE. The 
paper is organized as follows. In Sec. jlT we describe the 
NLSE used to study the UFG. In Sec~|lT| we study the 
dynamics of a single vortex in cigar shaped UFG. In Sec. 
|IV[ we study the CI in a pair of vortices which are gen- 
erated in the UFG by imprinting the appropriate phase. 
In the following Sec. [Vj we investigate the generation 



of vortex-antivortex pair by a Gaussian obstacle poten- 
tial, which on subsequent evolution decays via CI into 
vortex rings and sound. We conclude by presenting the 
summary of the results in the last section. 



II. NON-LINEAR SCHRODINGER EQUATION 
FOR UNITARY FERMI GAS 



The order parameter or the condensate wavefunction 
for the superfluid two component fermionic system, ex- 
cluding the normalization factor, is 



(1) 



where V't(i"j t) and ipiir, t) are the fermionic annihilation 
field operators for f and I spin states respectively. The 
order parameter can be interpreted as the wavefunction 
of the macroscopically occupied two particle state. The 
order parameter is normalized to total number of con- 
densed Cooper pairs [T9 ] . [45H47] . As mentioned earlier, 
the s-wave scattering length between two component spin 
states a — )■ — oo in UFG. The Lagrangian of the UFG 
with energy E and consisting of equal number of two 
spin states at T = K is [2l],[24], [^, [21], [30], [35] 



The energy of the UFG is given by 



E. 



(2) 



E = 



4^1^^'^ 
V{r,tM^ 



^^(37,2)2/3|^|10/3 

5to 



dr. 



(3) 



where m is the atomic mass of the fermionic species, ^ — 
0.44 [3^, and V{r, t) is the trapping potential. Using the 
action principle 



5 1 Ldt = 0, 
/ti 



(4) 



the time dependent NLSE describing the UFG at zero 
temperature is 



fi2 9fe2 



ih 



9V(r,i) 



dt 



(5) 



where mp — 2m is the mass of the fermionic pair. The 
order parameter ip satisfies the normalization condition 



\ip{r,t)\^dr = N. 



(6) 
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In the present work, we consider cigar-shaped trapping 
potential 



-(x +y + a z ), 



(7) 



where lo is the radial trapping frequency, a < 1 is the 
ratio of axial to radial trapping frequency. We consider 
UFG of "^K with N = 3500, w = lOOHz, and a = 0.2. In 
Refs. |33l|34], effective one dimensional (ID) NLSE for 
the cigar-shaped Fermi superfluids in the unitary limit 
has been proposed. The ID NLSE provides a good ap- 
proximation to the real three dimensional (3D) systems, 
provided the dynamics along radial direction is frozen. 
The CI instability being a three dimensional instabil- 
ity can not be studied using quasi-lD equation. Hence, 
in the present work, we consider the three dimensional 
equation and solve it numerically using the split time 
step Crank-Nicolson method [15]. In this method, the 
NLSE is discretized in space and time using the finite 
difference scheme. The time iteration is performed by 
splitting the Hamiltonian into two parts, one containing 
the spatial derivative part and the other containing the 
rest of the terms. The error involved in the splitting of 
the the Hamiltonian is proportional to the square of the 
time step. 

We can rewrite the Eq. ([s]) in scaled units using the 
transformations 



r = r'ansr,, t = t'oj \ 



^(r,<) 



3/2 
'lose 



(8) 



where the primed quantities are in scaled units and Oosc = 
^Jh/{mpU}) is the oscillator length. After dropping the 
primes, the scaled NLSE describing the superfluid UFG 
is 



V2 



2e(37r27V)2/3|0|4/3 + y(r,t) 



r,t) 

d(t>{Y,t) 
dt ' 



(9) 



This equation is used to study the dynamics of UFG in 
the present work. We use Aa; = Ay = Az = 0.15 and 
At = 0.001 as the spatial and time step sizes. The num- 
ber of time iterations used in the present work is 55000. 
In terms of unsealed units, the spatial and time step sizes 
are 0.169 /im and 1.59 x 10"'^ ms respectively. 



along z-axis. The stability of the vortex line in BECs 
can be inferred from the normal modes, which can be 
obtained by solving the equation obtained after lineariz- 
ing the Gross-Pitaevskii equation about the stationary 
configuration. In oblate condensates with a vortex line 
near z-axis there is one normal mode with negative fre- 
quency and it gives the precession frequency of the vor- 
tex in trap. On the other hand, the appearance of more 
than one normal modes with negative frequency in cigar 
shaped condensates, makes the vortex line increasingly 
susceptible to bending[5l]. In case of cigar-shaped UFG, 
a similar instability against bending can aid the reconnec- 
tion between the vortex and the antivortex line. Hence, 
before studying the dynamics of vortex-antivortex pair, 
it can be quite instructive to study the dynamics of a 
single vortex line in prolate UFG. 

We numerically generate the vortex line by phase- 
imprinting method, in which the order parameter is 
evolved in imaginary time under the constraint that it 
has the phase consistent with the presence of the vortex 
line at all times. To do so, while evolving the NLSE in 
imaginary time r = it, the order parameter (/)(r, r) is 
redefined as 



(r, T -I- 5t) — 0(r, r) exp 



tan 



-1 y-yo 



(10) 



after each iteration in imaginary time with the temporal 
step size of St. Here, (a;o, i;o) is the point of intersection of 
the vortex line with xy plane, and tan^^ [(y— j/o) /(s;— a^o)] 
is the azimuthal angle measured with the origin shifted 
to (xqjJ/OiO). This redefinition of the order parameter 
ensures that the solution of the NLSE always has vor- 
tex line passing through (a;o,2/o,0) and oriented parallel 
to z-axis. The solution obtained by this method is then 
evolved in real time to study the dynamics of the sys- 
tem. We first consider the dynamics of the UFG with 
an axial vortex line. In this case, the vortex line remains 
straight, and hence is stable against the bending insta- 
bility as is shown in Fig. [l] This is not the case for an 
off-axis vortex line, which tends to bend as it traverses 
across the condensate. This is evident from Fig.|2j where 
the upper and lower panel show the dynamics of the vor- 
tex line initially located at x = 1.0 Oosc and x = 3.0 Oosc 
respectively. Hence, like in BECs, an off-axis vortex line 
is unstable against bending in cigar-shaped UFG. 



III. DYNAMICS OF A SINGLE VORTEX LINE 



IV. CI IN VORTEX-ANTIVORTEX PAIR 
GENERATED BY PHASE-IMPRINTING 



In case of BECs, the dynamics of single vortex line in 
prolate condensates was studied experimentally in Ref. 
|49j . It was observed that vortex line in most cases 
was bent as was predicted in earlier theoretical studies 
[5DH55] . It was shown by Svidzinsky et al. |S1] that 
a straight vortex line is stationary only if it is oriented 



We can extend the phase- imprinting method, described 
in the previous section, to the vortex-antivortex pair. 
In order to numerically imprint a vortex-antivortex pair, 
phase contribution from both the vortex and antivortex 
has to be imprinted on the order parameter after each 
iteration in imaginary time. This can be ensured by re- 
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FIG. 1. The time evolution of the UFG with a single axial 
vortex line. The images are the isosurfaces corresponding to 
|<^(r)| = 0.012 aoBc^- Images with the labels (a) and (b) 
are the isosurfaces after ms and 79.6 ms of time evolution 
respectively. The origin of the coordinate system is located at 
the center of the box in each image. The dimensions of each 
box are 12.15 aosc x 8.55 Oosc x 54.15 aosc 



defining order parameter 0(r, t) as 



i tan 



(r, T + St) = (f>{r, r) exp 

X exp —i I tan 



-iV-yi 



X ~ Xi 

-1 y - 2/2 



(r, r) exp 



i tan 



X- X2 



X — Xi 



tan 



y-y2 

X - X2 



(11) 



after each iteration in imaginary time with the tempo- 
ral step size of St. Here tan~^[(y — yi)/{x — xi)] and 
tan~^[(y — y2)/{x — X2)] are the azimuthal angles mea- 
sured with the origin shifted to {xi,yi,Q) and (x2,?/2,0) 
respectively. The converged solution obtained by this 
method will have a vortex and antivortex at (xi, j/i) and 
{x2, 2/2) respectively. This solution is then evolved in real 
time to study the CI. As a case study, we imprint the 
vortex and antivortex at (— flosci —aosc) and (flosc, — Oosc) 
respectively. The stationary solution obtained by imag- 
inary time propagation is shown in Fig. [sjja); here the 
distance between the two parallel vortex lines is equal to 
2aosci ^'iid the origin of the coordinate system is located 
at the center of the box. The solution is then evolved 
in the real time using Eq. (|9]). During the initial stages 
of the evolution vortex-antivortex follow the trajectories 
which is reminiscent of the trajectory followed by a vor- 
tex dipole in the pancake-shaped traps [H]. To be pre- 
cise, the vortex moves in the region to the left of x = 
plane, while the antivortex moves in the region to the 
right of a; = plane. In order to understand it, the 
vortex-antivortex pair at i = ms can be considered to 




FIG. 2. Upper panel shows the dynamics of the UFG with an 
off-axis vortex line initially located at a; = 1.0 Ooac, whereas 
the lower panel shows the dynamics with an off-axis vortex 
line initially located a,t x — 3.0 flosc. The images are the 
isosurfaces corresponding to |0(r)| — 0.012 aosJ"^ ■ In upper 
panel, images with the labels (a) and (e) are the isosurfaces 
after ms, 23.8 ms, 47.7 ms, 63.7 ms, and 79.6 ms of time 
evolution respectively. In the lower panel, images with the 
labels (f) and (j) are the isosurfaces after ms, 12.7 ms, 47.7 
ms, 63.7 ms, and 82.8 ms respectively. The origin of the 
coordinate system is located at the center of the box in each 
image. The dimensions of each box are 12.15 floac x 8.55 Oosc x 
54.15 aosc 



be equivalent to a highly stretched vortex ring with cen- 
ter at (0, —aosc)- The motion of this vortex ring will be 
due to the contribution from two velocity fields: a) the 
self induced velocity which will propel the ring in pos- 
itive y direction and b) the precession induced by the 
inhomogeneity of the condensate. Initially the self in- 
duced velocity dominates and the ring travels forward, 
later on the precession effects dominate causing the ring 
to move backward. The precession will result in the ex- 
pansion of the ring along radial direction for y > and 
contraction for j/ < 0. After the ring returns close to 
its initial position, each element of the ring (vortex and 
antivortex) has developed a curvature due to aforemen- 
tioned expansion as is shown in Fig. [sj^b). As the pair 
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FIG. 3. (Color online) The time evolution of the cigar-shaped 
UFG with vortex-antivortex imprinted in it at t = ms. The 
images are isosurfaces corresponding to |(;/)(r)| — 0.012 ao^J^ ■ 
In the upper panel, images with labels (a) to (e) are the iso- 
surfaces after ms, 20.7 ms, 28.6 ms, 31.8 ms, and 35.0 ms of 
time evolution respectively. In the lower panel, images with 
labels (f) to (j) are the isosurfaces after 38.2 ms, 41.4 ms, 
47.7 ms, 54.1 ms, and 79.6 ms of time evolution respectively. 
The origin of the coordinate system is located at the center 
of the box in each image. The dimensions of each box are 
12.15 Uosc X 8.55 Uosc X 54.15 aosc. 



moves along positive y direction, the distance between 
the vortices at their edges decreases further causing the 
first reconnection event. This results in an actual vortex 
ring as is shown in Fig. |3][c) . The first reconnection event 
is aided by bending instability of the off-center vortices. 
Due to this bending instability, vortex and antivortex 
bend in opposite direction resulting in the decreased sep- 
aration between them at the edges [see Fig. [Sjjb)]. The 
reconnection event is also accompanied by the emission 
of the sound as is evident from ripples on the isosurface 
in Fig. [sj^c). Also, in the presence of trapping poten- 
tial, parts of vortex ring in higher density region move 
slower than those in lower density regions [M] . This leads 
to the differential velocity field experienced by the ring. 
The differential velocity field causes the ring to bend at 
the ends [Fig. [sj^c)]; this in turn generates the Kelvin 
waves which aid the growth of Cf as is shown in Figs. 



^c) and (d). For a sufficiently large amplitude of these 
Kelvin waves, the second reconnection event takes place 
at 31.8 ms [Fig. [Sjjd)]. This results in the generation of 
two daughter vortices and one parent vortex [Fig. [sje)]. 
Reconnection events are typically accompanied by the 
generation sound waves and further generation of Kelvin 
waves. With further evolution, the parent ring under- 
goes the same bending at the edges [Fig. [sjf)], followed 
by the generation of second set of daughter rings [Fig. 
|3][g)]. As more time elapses, more reconnection events 
and the penetration of the vortex rings in the UFG's 
surface take place [Figs. [Sf^h-j)]. At the end of the time 
evolution, the initial vortex-antivortex pair has mostly 
decayed into sound waves due to CI as is shown in Fig. 
ij)- 



V. CI IN VORTEX-ANTIVORTEX PAIR 
GENERATED BY A MOVING OBSTACLE 
POTENTIAL 

Experimentally vortex-antivortex pairs have been cre- 
ated in the BEC by moving the condensate across a re- 
pulsive Gaussian obstacle [15]. The generation of vor- 
tex rings by moving a spherical obstacle across a UFG 
has also been studied theoretically |20]- The method of 
generating vortex antivortex pairs by a moving repulsive 
potential is also a good prospect to experimentally ob- 
serve CI in BECs [12]. Taking all these into account, 
we also study the evolution of UFG under the influence 
of a moving obstacle potential. We employ an obstacle 
potential 

K..(t)^V-oW^exp|-2 -^ + [^-^f(^)n , (12) 

created by Gaussian laser beam. Here wq is the waist size, 
Wz = wo\/r^-^V^; where is the Rayleigh range, 
is beam width or spot size, and Vo{t) is the maximum 
potential created by the beam at the center of beam 
waist which lies at (0,?/o)- We use the same obstacle 
potential parameters as used in Ref. [T^] : Vq — 25 huj, 
Zfl — 22 /im, and wq = 2.2 iim. The stationary state iso- 
surface with the obstacle potential located at (0, — Oosc) 
is shown in Fig. |4][a); the figure clearly shows the de- 
pletion of UFG in the region of the obstacle potential. 
We move the obstacle with the speed of 500 /im/s along 
y-axis. This speed is a fraction of the speed of sound in 
UFG at the trap center. We estimate the speed of sound 
at the trap center by using the expression for speed of 
sound in a homogeneous UFG, i.e., c = y/^vp/VS, where 
vp = 2(37r^A^|(/)|^)^/^ is the Fermi velocity for ideal Fermi 
gas [TS]. Using this expression, the speed of the sound 
at the trap center is ~ 1800 ^m/s. While moving the 
obstacle its amplitude is linearly ramped down. The rate 
of ramping down is such that the obstacle potential be- 
comes zero at (0,4 Oosc)- 

As the obstacle moves along positive y direction, it ere- 
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FIG. 4. (Color online) Evolution of the UFG when the ob- 
stacle initially located at (0, — Oosc) is moved with the speed 
of 550 nm/s. The images are the isosurfaces corresponding 

to \4>{t)\ = 0.012 aosJ"^ ■ In the upper panel, images with the 
labels (a) to (e) are the isosurfaces after ms, 6.4 ms, 20.7 
ms, 25.5 ms, and 28.6 ms of time evolution respectively. In 
the lower panel, images with the labels (f) to (j) are the iso- 
surfaces after 30.2 ms, 35.0 ms, 38.2 ms, 43.0 ms, and 82.8 ms 
of time evolution respectively. The origin of the coordinate 
system is located at the center of the box in each image. The 
dimensions of each box are 12.15 Oosc x 8.55 ttosc x 54.15 flosc- 



tion events [Fig. |4][e)] With further time evolution, CI 
leads to more reconnection events, penetration of the vor- 
tices in the UFG's surface, and generation of sound waves 
[Figs. |4f^f-j)] as happened in the case of phase- imprinted 
vortex-antivortex pair. At the end of the evolution pe- 
riod, only sound waves are discernible in the system [Fig. 

iJ)]- 



VI. SUMMARY AND CONCLUSIONS 

We have numerically investigated the CI in the cigar- 
shaped superfluid Fermi gas in the unitary limit. We 
use NLSE based on the time dependent density func- 
tional theory to study CI. We numerically solve this 
equation using split time step Crank-Nicolson method 
to study the dynamics UFG with a single vortex line 
as well as a vortex-antivortex pair. In order to gener- 
ate the vortex-antivortex pair, we employ two methods: 
(a) phase-imprinting and (b) moving a Gaussian obsta- 
cle potential across the system. We have shown that a 
single vortex line in cigar shaped UFG is stable only if it 
is oriented along axial direction. On the other hand, an 
off-axis vortex line tends to bend as it traverses across 
the condensate. We have shown that a vortex-antivortex 
pair oriented along the axial direction in a cigar-shaped 
UFG is susceptible to the CI mechanism. We observe 
that frequent vortex reconnections, inhomogeneity gen- 
erated Kelvin waves, and sound generation characterize 
the CI in an inhomogeneous UFG like in a BEC. We 
have shown, using the experimentally realizable param- 
eters for trapping potential and number of atoms, that 
it should be possible to observe the signatures of CI in 
a vortex-antivortex pair generated by moving a blue de- 
tuned laser beam across the UFG. 
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